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We consider quantum charge pumping of electrons across a superconducting double barrier struc- 
ture in the adiabatic hmit. The superconducting barriers are assumed to be reflect ion- less so that 
an incident electron on the barrier can either tunnel through it or Andreev reflect from it. In this 
structure, quantum charge pumping can be achieved (a) by modulating the amplitudes, Ai and A2, 
of the gaps associated with the two superconductors or alternatively, (b) by a periodic modulation 
of the order parameter phases, (pi and (j)2 of the superconducting barriers. In the former case, we 
show that the superconducting gap gives rise to a very sharp resonance in the transmission resulting 
in quantization of pumped charge, when the pumping contour encloses the resonance. On the other 
hand, we find that quantization is hard to achieve in the latter case. We show that inclusion of 
weak electron-electron interaction in the quantum wire leads to renormalisation group evolution of 
the transmission amplitude towards the perfectly transmitting limit due to interplay of electron- 
electron interaction and proximity effects in the wire. Hence as we approach the zero temperature 
limit, due to renormalisation group flow of transmission amplitude we get destruction of quantized 
pumped charge. This is in sharp contrast to the case of charge pumping in a double barrier through 
a Luttinger liquid where quantized charge pumping is actually achieved in the zero temperature 
limit. 

PACS numbers: 73.23.-b,74.78.Na,74.45.+c 



I. INTRODUCTION 

The phenomena of quantum charge pumping corre- 
sponds to a net flow of DC current between different 
electron reservoirs (at zero bias) connected via a quantum 
system whose system parameters are periodically mod- 
ulated in time^i^. The zero bias current is obtained in 
response to the time variation of the parameters of the 
quantum system which explicitly break the time rever- 
sal symmetry. It is worth mentioning that breaking of 
the time-reversal symmetry is necessary but not a suf- 
ficient condition in order to get net pumped charge in 
unit cycle. For obtaining a net pumped charge, par- 
ity or spatial symmetry must also be broken. Within a 
scattering approach, if the time period of modulation of 
the scattering system parameters is much larger than the 
time the particle spends inside the scattering region, adi- 
abatic limit is reached. In this limit, the pumped charge 
in a unit cycle becomes independent of the pumping fre- 
quency. This is referred to as "adiabatic charge pump- 
ing" . Experimentally charge pumping has been observed 
in mesoscopic systems involving quantum dots and car- 
bon nanotubes"^'^'^ . Ofcourse one has to be very careful in 
interpreting the experimentally observed pumped charge 
as it can be faked by rectification effects as was pointed 
out by Brouwer— . 

In the recent years, there has been an upsurge of re- 
search interest in exploring the effects due to inclusion 
of electron-electron (e — e) interaction on the pumped 



charge^i^i^iiSiiiiii^. In this article, we explore the effect 
of inter-electron interaction on the charge pumped across 
a superconducting double barrier (SDB) system^^ in the 
context of one-dimensional (1-D) quantum wire (QW). 
Pumping of free electrons across 1-D quantum well was 
studied earlier in Refs. Hflalla . where using Brouwer's 
formula^, it was shown that the pumped charge can be 
expressed as a sum of two contributions, viz., a dissi- 
pative part and a quantized topological part, the latter 
being independent of the details of the pumping con- 
^Q^piTjis rpj^g dissipative part was found to be pro- 
portional to the conductance through the system on the 
pumping contour in the parameter space while the topo- 
logical part was non-zero only if the pumping contour 
enclosed a resonance. Hence in order to obtain quan- 
tized pumped charge, one needs to reduce the dissipa- 
tive part as much as possible. This is very naturally 
achieved if one considers pumping through a quantum 
well in a 1~D interacting electron gas^ii^ (Luttinger liq- 
uid) as in this case interaction correlations make the res- 
onance very sharp thereby reducing the conductance on 
the contour enclosing the resonance to vanishingly small 
values in the zero temperature limit. This leaves be- 
hind a quantized topological part. The pumped charge 
was shown to converge to a quantized value asymptoti- 
cally. This was obtained using a perturbative approach 
for the case of a weakly interacting electron gas followed 
by "Poor-man's scaling approach"^". In this article, we 
show that the presence of inter-electron interaction in 



the wire leads to nontrivial scattering processes due to 
proximity effects wliicli leads to a power law reduction of 
the pumped charge from the quantized value (as opposed 
to enhancement) in the adiabatic limit as one lowers the 
temperature. Quantum charge pumping using various 
setups involving superconductor has been a topic of ma- 
jor interest in recent pas t^-'^'^^i^^'^"^i^^'^^i^''''^^i^^i'^'-' . Specif- 
ically we consider pumping of electrons (in the adiabatic 
limit) across a SDB structure, as depicted in Fig.[T] Till 
date no experiment has been carried out in the context 
of charge pumping for the case of superconducting bar- 
rier. Experimentally it might be possible to design a 
SDB structure by depositing thin strips of superconduct- 
ing material on top of a single ballistic QW (like carbon 
nanotubes) at two places, which can induce a finite su- 
perconducting gap in the barrier regions of the QW as 
a result of proximity of the superconducting strips. In 
our simple-minded theoretical modelling of the system 
we assume that the superconducting barrier (SB) to be 
reflection-less so that an incident electron on the barrier 
can either tunnel through it or Andreev reflect from it. 
Within the simplified theoretical model, we explore two 
scenarios to achieve quantization of pumped charge — (a) 
by periodic modulation of amplitudes Ai and A2 of the 
gap at the two SB or alternatively, (b) by periodic modu- 
lation of the order parameter phases (f>i and (f>2 associated 
with the two SB. For free electrons in the QW we show 
that in the Ai — A2 plane, there is an isolated sharp 
resonance point in transmission probability across the 
SDB structure. On the other hand transmission probabil- 
ity across the double barrier has a line of sharp resonances 
in the (j)i —<t>2 plane. As mentioned earlier, in order to ob- 
tain quantized pumped charge, the transmission on the 
pumping contour should be as small as possible. When 
we consider Ai and A2 as the pumping parameters, we 
can always choose a pumping contour which completely 
encloses the isolated resonance and hence it is possible to 
achieve quantization of charge if the resonance is sharp 
enough. However, in the 0i — 02 plane, we have a line of 
resonances. Any closed contour enclosing the resonances 
will surely cross the resonance line at least twice thereby 
increasing the dissipative part and consequently resulting 
in destruction of quantization of pumped charge. Inter- 
estingly enough, inclusion of weak e — e interaction in the 
wire results in a RG flow of the transmission amplitude 
towards perfectly transmitting limit due to proximity in- 
duced effect on the interacting electrons in the QW as we 
lower the energy scale such as temperature. Hence the 
sharpness of resonance is lost due to RG enhancement of 
transmission through individual SB resulting in complete 
destruction of quantized charge pumping as we go down 
in temperature. It is worth noticing that the consequence 
of inclusion of correlations due to e — e interaction is just 
opposite here with respect to the case of double barrier 
in a Luttinger liquid^. 

This article is organised as follows. In Sec. II, we discuss 
the modelling of SDB in a 1-D QW and calculate the 
transmission and Andreev reflection (AR) amplitudes of 
the system. In Sec. HI, we discuss the renormalisation 
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FIG. 1: A one-dimensional quantum wire (e.g. carbon nan- 
otube) connected to two reservoirs, labelled by 1 and 2. The 
two thin strips on the wire depict layers of superconducting 
material deposited on top of the wire. The superconducting 
strips are connected to contacts labelled 3 and 4. 



group (RG) flow for transmisson and AR for the SB. In 
Sec. IV, we discuss our RG scheme for the SDB and cal- 
culate the pumped charge. In the end, we discuss our 
results and give conclusions in Sec. V. 



II. SUPERCONDUCTING BARRIER 

Quantum transport in SB structure was considered in 
past in Ref. il3i . Here we consider a very similar set-up 
comprising of a ballistic 1-D QW with two short, but 
finite superconducting patches as shown in Fig. [1] Here, 
A'*^ is the pair potential on the two patches {i refers to 
the index of the strip). Following Ref. [ij, we use the 
Bogolubov— de Gennes (BdG) equation^i^^ to calculate 

the transmission amplitude lie and the AR amplitude, 
Tg^ , where i is the barrier index. The space dependance 
of the order parameter (which also acts as the scattering 
potential) for the incident electron can be expressed as 

V{x) = AWe''^ie(a;)e(-x + a) + A«e*'^^ 

Q[x-{a + L)]Q[-x+{2a + L)] (1) 

where, a is the width of the SB and L is the distance 

between the two barriers. 

Hence the BdG equations can be written as, 
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Solving the BdG equation in the normal and supercon- 
ducting regions and matching the solution at x = and 
X = a, we get 
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= ^2to(£;f±(£;2-a2W)i/2), m^ = 

,(A(*)/£;)2)i/2)]i/2 Here, 41 1«, r«, ri^^ are 

the transmission and AR ampHtudes. m is the effective 
mass of the electron in the wire, Ep is the Fermi energy 
for the electrons in the superconducting region, and E is 
the Fermi energy of electrons in the normal region of the 
QW, measured with respect to Ep. Hence the scattering 
matrix for the single SB problem for an incident electron 
or hole is given by 
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Using the 5- matrix given by Eq.[5l we obtain the effective 
S'-matrix for the double barrier system^. We assume 
particle- hole symmetry, hence t^e = thh and r^h = fhe- 
The net transmission and net AR amplitude through the 
double barrier are 
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where hq = ^1ra{Ep ± _B). In order to obtain quan- 
tization of pumped charge, we choose to operate in the 
sub-gap regime i.e., E < A. In this regime, |TeeP has 
sharp resonances at discrete values of E/A for a given 
value of (pi — (j)2^^-- These resonances result from multi- 
ple AR of electron to hole and vice- versa inside the double 
barrier. 



III. 



WIRG STUDY OF JUNCTIONS 



We study the effects of inter-electron interactions in 
the wire on the S'-matrix characterizing the supercon- 
ducting barrier using the RG method introduced in 
Ref. 1 201. and the generalizations to multiple wires in 
Refs. [34II35I The basic idea of the method is as fol- 
lows. The presence of back-scattering (reflection) induces 
Friedel oscillations in the density of non-interacting elec- 
trons. Within a mean field picture for weakly interacting 
electron gas, the electron not only scatters off the po- 
tential barrier but also scatters off these density oscilla- 
tions with an amplitude proportional to the interaction 
strength. Hence by calculating the total reflection ampli- 
tude due to scattering from the scalar scatterer and from 
the Friedel oscillations created by the scatterer, we can 
include the effect of e — e interaction in calculating trans- 
port. This can now be generalized in a similar spirit to 



the case where there is, besides non-zero reflection also 
non-zero AR which turns an incoming electron into an 
outgoing hole due to proximity effects as done in Ref. [Sa 
and then generalized to multiple wire superconducting 
junction in Ref. 37. 

The fermion field on each wire can be written as. 
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where i is the wire index, s is the spin index which can 
be 1,1 and /, O stands for outgoing or incoming fields. 
Note that "^ i{x){'^ oi^)) are slowly varying fields on the 
scale of kp and contain the annihilation operators as well 
as the slowly varying wave- functions. The expectation 
values for the density ( ^]^j,^is ) gives (dropping the wire 
index) , 
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and (■0l|^OT> = (V^lj,'0Oi) 
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Hence, besides the density, the expectation values for 
the pair amplitude ( ^Is^lg ) and its complex conjugate 
( ^is^is ) are also non-zero and are given by (dropping 
the wire index) 
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So, we see that the Boguliobov amplitudes fall off as l/x 
just like the normal density amplitudes. 
We now allow for short-range density-density interactions 
between the fermions 
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to obtain the standard four-fermion interaction Hamilto- 
nian for spin-full fermions as 

Hint = Jdx[gi(^^\^^l^.-^j^^oi + *n*oi*/i*oi 
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(13) 

where gi and ^2 are the interaction parameters^. 

The effective Hamiltonian can be derived using a 
Hartree— Fock (HF) decomposition of the interaction 
Hamiltonian. The charge conserving HF decomposition 



leads to the interaction Hamiltonian (normal) of the fol- 
lowing form on each half wire, 
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(We have assumed spin-symmetry i.e. r^ — r^ — r.) 
This has been derived earlier—. Using the same method, 
but now also allowing for a charge non-conserving H F de- 
composition we get the (Andreev) Hamiltonian 
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The e — e interaction induced amplitude to go from 
an incoming electron wave to an outgoing electron wave 
under e"*^'"** (for electrons with spin) is given by^ 
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where Oi= {gi — '2,gi)l2-KhvF and d is the short distance 
cut-off for the RG flow. Analogously, the amplitude to 
go from an incoming electron wave to an outgoing hole 
wave under e^^'^'-"t* is given by^ 
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where a' = {gi + g2)/2'!rhvp. 

These logarithmic corrections to the bare reflection am- 
plitude and the AR amplitude can be summed up using 
a "Poor-man's scaling approach" ^^ which flnally leads a 
RG equation for r and r^. 



IV. RENORMALISATION GROUP SCHEME 
AND THE PUMPING FORMULA 

We include the effects due to proximity of supercon- 
ductor and e — e interaction in the wire via a RG ap- 
proach developed very recently^I for the case of 1-D nor- 
mal metal— superconductor— normal metal (NSN) junc- 
tion. As we are only interested in the coherent regime 
{Lt >> L, Lt is the thermal length), we can effectively 
treat the SDB system (NSNSN junction) as a single bar- 
rier (NSN junction) as far as RG is concerned. 

Hence the effective two-channel S'-matrix for this dou- 
ble barrier system can be written as 
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where all the amplitudes and phases associated with the 
matrix elements are functions of the time- varying param- 
eters, Vi{t) = Vb -I- Pcos(cjt + {-ly^'^ri) where i = 1, 2 



stands for the barrier index. Vi = A^ and Vi = (pi are the 
two possible pumping parameters. The reflection coeffi- 
cients are not the same (phases can differ) because the 
time-varying potentials explicitly violate parity. In prin- 
ciple the S'-matrix also violates time-reversal invariance. 
But since in the adiabatic approximation, we are only 
interested in instantaneous hamiltonian. Note that the 
instantaneous S'-matrix can mimic a time-reversal sym- 
metric S'-matrix. 

Using the modified Brouwer's formula^, the pumped 
charge can directly be obtained from the parametric 
derivatives of the S-matrix elements. It is worth men- 
tioning that even though Brouwer's formula is valid for 
non-interacting electron system, we are able to use it here 
because effects due to interactions in the wires can be 
taken care of by the renormalization of the bare S'-matrix 
obtained for the free-electron case. 

For single channel S-matrix, we have 
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where Sij denote the elements of the S'-matrix. 



(19) 



Note 



the negative sign in the above expression, which results 
from the fact that S'n corresponds to conversion of an 
electron into a hole. Thus, the pumped charge is directly 
related to the amplitudes and phases that appear in the 
S-matrix. Inserting Eq. [18] in Eq. [191 
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Here G{t) — |Tee(i)|" is the instantaneous two terminal 
linear conductance (labelled by 1,2 in Fig. [T|), in units 
of 2e2//i. The first term on the RHS in Eq. ^ is clearly 
quantized since e*^ returns to itself at the end of one cy- 
cle. So the only possible change in in a period can 
be in integral multiples of 2tt i.e., 9{t) — > 9{Q) -t- 27rn 
[n -^ integer). The second term is the 'dissipative' 
term which prevents the perfect quantization. The sec- 
ond term is directly proportional to the two terminal 
Landauer— Buttiker conductance for the system on the 
pumping contour. The relative sign between 9 and in 
the expression for pumped charge in Eq. [201 originates 
from the AR process, which converts an electron to a 
hole. This is in contrast to what has been found for 
the normal double barrier problem^^. For a reflection- 
less junction, the basic idea of the RG method is as fol- 
lows. The presence of a superconductor induces a fl- 
nite yet weak pair potential in the QW resulting in scat- 
tering of incoming electrons to outgoing holes (Andreev 
processes) in the wire, away from the junction. Hence 
by calculating the total AR amplitude, due to scattering 
from the NSN junction and the (weak) pair potential in 
the wire perturbatively in interaction strength and fol- 
lowed by "Poor-man's scaling" approach, we obtain the 
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FIG. 2; Contours of transmission probability, \Tee\ in the Ai — A2 plane at two different values of length scale, Lp = 1 and 
Lp = 10, at which the RG flow is cut-off for values of 1^(0) = 0.1 and V{2kF) = 0.1. The red ellipse Ci represents the pumping 
contour. 
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FIG. 3: Contours of transmission probability, |Tee| in the 01 — (^2 plane at two different values of length scale, Lp — 1 and 
Lp = 10 at which the RG flow is cut-off for values of 1^(0) = 0.1 and V{2kF) — 0.1. The red circle C2 represents the pumping 
contour. 



RG equation for the elements of the effective S'-matrix of 
the SDB structure in the coherent regime (Lr > L). So, 
the entries of 5'-matrix therefore become functions of the 
length scale Lp due to the RG flow. The RG flow can 
also be considered to be a flow in the temperature since 
the length scale Lp can be converted to a temperature 
scale using the thermal length Lt — hvp/ikBT). Hence, 
the RG flow has to be cut-off by either Lt, or the system 
size Ls, whichever is smaller—. 

Without loss of generality, we can calculate the renor- 
malized 5'-matrix at different length scales or equiv- 
alently at different temperatures at any point on the 



pumping contour. Hence, to avoid unnecessary compli- 
cations arising due to the RG flow of phases associated 
with S'-matrix elements {6,d',(j>), we choose to calculate 
the RG flow of the S'-matrix when the barriers are sym- 
metric. This symmetry leads to vanishing of the RG flow 
of the phases hence making the calculation algebrically 
simple. 

The RG flow of the normal transmission (and AR) am- 
plitudes and phases are^ 
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FIG. 4: Pumped charge Q, for pumping in Ai — A2 plane, is 
shown as a function of the dimensionless parameter I where 
I — ln{Lp/d) and Lp is either Lt = hvp/ksT at zero bias or 
Lv = hvp/eV at zero temperature and d is the short distance 
cut-off for the RG flow. The three curves correspond to three 
difi'erent values of V^(O) and V{2kp). 
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Here I — ln{Lp/d) where d is the short distance cut-off 
for the RG flow and we have considered the fully symmet- 
ric case, i.e. 9 = 0'. Unitarity of the S'-matrix in Eq.fTSl 
implies that (p — 9 = tt/2 + 2mT (n -^ integer). This 
simplifies the equations for RG flow for the AR amplitude 
and phase, 
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Here, a' = (52 +5i)/27r?iwi^ where 51,52 are the running 
coupling constants whose bare values are set by 51 {Lp = 
d) = V{2kF) and g2{Lp = d) = V{0); V{x) being the 
inter-electron interaction potential. We now integrate 
the RG equation for T^e complimented by the RG flow of 
51 and 52^^ to obtain the Lp dependence of T^e as 



(23) 



Here T^^ and R'^i^ are the values of T^e and Reh at length- 
scale L and ai = V{0)/2TThvp and 02 = V{2kp)/2nhvp. 
There are two points worth mentioning here : (a) the 
transmission increases with increasing Lp which is a con- 
sequence of the fact that the proximity effect due to su- 
perconductor induces an effective attractive interaction 
between the electrons, hence rendering the (Andreev) 
back-scattering an irrelevant operator, and (h) the ex- 
pression for TeeXLp) is not in the form of a pure power 
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law even at T^^ -^ limit, as is expected from Luttinger 
Liquid physics because of the RG flow of the g\ , 52 pa- 
rameters. Also, it is important to note that we take the 
short-distance cut-off d to be the distance between the 
two barriers (L) since this is the length scale at which 
we glued the two barriers to a single barrier as far as 
RG is concerned. Using this, we can obtain the scaling 
behavior of the pumped charge (Q) as a function of the 
length scale Lp (or the temperature T) . In terms of the 
Landauer— Buttiker conductance, Gq = (2e^//i) |T|?gp, 
using Eq. [231 we obtain the pumped charge as 
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'-9 + 4' and as earlier, Gq is expressed in unit of 
Qint is the integer contribution of the flrst term 
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FIG. 5: The plot shows the variation of the AR phase </> with 
time t, along the pumping contour Ci in the plane of Ai — 
A2 and the inset shows the variation of the same along the 
pumping contour C2 in the plane of <^i — 4>2. 



V. RESULTS AND DISCUSSIONS 

Pumping in the Ai — A2 plane : Here the pumped 
charge is obtained by periodically varying the top 
gate voltage which controls the Fermi energy of the 
electrons in the superconducting region. Hence it 
amounts to varying E/A for the two barriers pe- 
riodically. Just like the double barrier problem, in 
this case too we observe resonant transmission of 
electrons at discrete values of E/A for flxed val- 
ues of 4>i and 02 . These discrete values correspond 



to the existence of quasi-bound states formed in- 
side the SDB which are quite different from their 
normal double barrier counterpart as they are pro- 
duced due to superposition of both electron and 
hole states and not just any one of them. In Fig. [2] 
(left panel), we see sharp resonance in transmis- 
sion probability (|Teep) in the Ai — A2 plane for 
L = 1. We employ the solutions to the RG equa- 
tions (Eq. [2T|l to obtain the renormalized surface 
of transmission in the plane of Ai — A2 for a value 
of L = 10, this is shown in Fig. [2] (right panel). 
Note that the RG flow is such that the transmis- 
sion increases in the entire Ai — A2 plane, hence 
reducing the sharpness of resonance and resulting 
in an increase of transmission (conductance) on 
the pumping contour Ci giving rise to reduction 
in the pumped charge from its quantized value (see 
Fig. 2]). From Fig.[5l we notice that the AR phase (j) 
shows a total drop in its value by a factor of 27r dur- 
ing its time evolution along the contour Ci. This 
drop corresponds to the quantization of the topo- 
logical part in the expression for pumped charge Q 
(Eq. [2DD to the value of e. 
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FIG. 6: Pumped charge Q, for pumping in (jii — (^2 plane is 
shown in the figure as a function of the dimensionless param- 
eter I where I = ln{Lp/d) and Lp is either Lt = hvp /kpT 
at zero bias or Lv = hvp/eV at zero temperature and d is 
the short distance cut-off for the RG flow. The three curves 
correspond to three different values of V{Q) and V(2kF)- 



2. Pumping in the 0i — (/>2 plane : In contrast to 
the previous case, here we obtain two sharp lines 
of resonances for the transmission function in the 



4)1 — 4>2 plane. Again we observe in Fig. [3] that the 
RG flow (Eq. [2T|l results in reduction of the sharp- 
ness of the resonance. We consider a pumping con- 
tour C2 which encloses parts of both the resonance 
lines in the 0i — (/)2 plane. The intersection of the 
pumping contour C2 with the lines of resonance 
results in vanishing of the topological part. This 
can be seen by observing the time-evolution of the 
AR phase along contour C2 as shown in the inset 
of Fig. \5\ In this case the drops are exactly com- 
pensated by corresponding rises in phase (j) by same 
amount, leading to a net zero topological contribu- 
tion to the pumped charge. Hence for small values 
of Lp (see Fig. [5]), the pumped charge is almost 
zero. This is because the topological part is iden- 
tically zero while the dissipative part is non-zero 
but vanishingly small (due to the resonance being 
very sharp) as the conductance on most part of 
the contour is negligible. As we go to the larger 
Lp values, the pumped charge shows an interest- 
ing non-monotonic behavior, purely coming due to 
the variation of the dissipative part. 

In conclusion, we show that pumping in the Ai — A2 
plane is much more efficient as opposed to that in 0i — (j)2 
plane. We also demonstrate that the quantization of the 
pumped charge is lost in Ai — A2 plane if we include 
correlations due to proximity effects in the 1-D QW. Al- 
though if the barriers are reflecting then according to 
RG, the system will flow to the disconnected fixed point 
(r — 1) at low temperature. In that case the sharp trans- 
mission resonance would appear in the parameter plane 
of back-scattering strength of the first and the second 
barriers. If the pumping contour encloses the transmis- 
sion resonance, then in the zero temperature limit, the 
dissipative part of the pumped charge will become van- 
ishingly small resulting in quantized pumped charge. So 
for the SDB system with small normal reflection, pumped 
charge will eventually converge to a quantized value in 
the zero temperature limit. 
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